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Abstract 

This paper presents a new approach to the Hamiltonian structure of isomon- 
odromic deformations of a matrix system of ODE's on a torus. An isomonodromic 
analogue of the SU(2) Calogero-Gaudin system is used for a case study of this 
approach. A clue of this approach is a mapping to a finite number of points on 
the spectral curve of the isomonodromic Lax equation. The coordinates of these 
moving points give a new set of Darboux coordinates called the spectral Darboux 
coordinates. The system of isomonodromic deformations is thereby converted to 
a non-autonomous Hamiltonian system in the spectral Darboux coordinates. The 
Hamiltonians turn out to resemble those of a previously known isomonodromic sys- 
tem of a second order scalar ODE. The two isomonodromic systems are shown to 
be linked by a simple relation. 
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Comments on revised version 



This is a revised version of the paper published in Journal of Mathematical Physics vol. 
44, no. 2 (2003), pp. 3979-3999. It turned out after publication that the published 
version contains serious errors on the construction of the Lax and zero-curvature equa- 
tions. Namely, the M-matrices Mj(z) have to be corrected by an extra diagonal matrix 
diag(pj, — pj); the zero-curvature equations without these terms lead to a contradiction. 
The emergence of these correction terms was already pointed out in the paper of Ko- 
rotkin and Samtleben [10J. As they stressed therein, this problem stems from the fact 
that the isomonodromic system in question is a constrained system. The quantities pj, 
j = 1, . . . , N, may be interpreted as the Lagrange multipliers in the Hamiltonian formu- 
lation of this constrained system. This issue is also partly related to the treatment of the 
coefficient k that arises in the construction of spectral Darboux coordinates. 

To correct these errors, Section III has been fully revised. Firstly, the "fake" Hamil- 
tonians Tij are replaced by the Hamiltonians Tij with a correction term proportional to 
the left hand side of the constraint. The equations of motion of the Calogero variables 
q,p and the spin variables are defined by these Hamiltonians. The i-dependence of 
p/s is not determined at this stage. Secondly, the M-matrices Mj(z) are corrected by the 
diagonal matrix diag(pj, — pj). The Lax and zero-curvature equations are reformulated 
in terms of the corrected M-matrices Mj(z). The diagonal part of the zero- curvature 
equations then yields a set of new differential equations for pj's. As it turns out, the 
equations of motion of q,p and A^s ensure the integrability, in the sense of Frobenius, of 
these equations. Thus, as Korotkin and Samtleben [TO] noted, the t-dependence of p/s 
are eventually determined by the zero-curvature equations. 

The subsequent sections are left mostly intact, except for the last part of Section VI 
that contained a wrong statement on k. 

I would like to express my gratitude to Henning Samtleben for helpful suggestions on 
this issue. 
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I Introduction 



The notion of isomonodromic deformations originates in the celebrated work of R. Fuchs 
PP. Fuchs studied isomonodromic deformations of a second order linear ODE of the form 

d 2 y fab c d 3 

dz 2 + {~? + (z-1) 2 + {z - t) 2 + z(z - 1) ~~ A{z - A) 2 

tit^. | A(A - l)u \ 
z{z-l){z-t) z{z - \){z - X) ) y 

with five regular singular points z — 0, 1, oo, t, A on the Riemann sphere, and discovered 
a nonlinear ODE that is nowadays called the sixth Painleve equation. His work was 
soon generalized by Gamier [2] and Schlesinger [3] in two different directions. Whereas 
Gamier extended the work of Fuchs to a second order linear ODE with more singularities 
(including irregular ones as well), Schlesinger studied a matrix system of the form 

dz piz-tj 
and obtained the so called Schlesinger system 



dtj tj j t\ 

that characterizes isomonodromic deformations. It turned out afterwards [3] that Gar- 
nier's isomonodromic deformations with an arbitrary number of regular singular points 
can be reconstructed from the 2x2 Schlesinger system. 

The next stage of generalization is, naturally, isomonodromic deformations on a torus. 
This issue was first tackled by Okamoto [5 , who obtained a system of isomonodromic 
deformations of a second order scalar ODE on a torus. One of his remarkable results is 
that the isomonodromic system on a torus can be formulated as a Hamiltonian system 
in much the same way as Garnier's isomonodromic system on a sphere was converted to 
a Hamiltonian system Iwasaki [7j extended Okamoto's work to scalar ODE's of an 
arbitrary order on an arbitrary compact Riemann surface, and elucidated the geometric 
origin of the Hamiltonian structure that Okamoto derived. The study of isomonodromic 
systems on a torus was further refined by Okamoto himself [Hj and Kawai 0. 

As regards matrix systems, Korotkin and Samtleben [TD] constructed an example of 
isomonodromic deformations of a 2 x 2 matrix system on a torus. Levin and Olshanetsky 
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[TT] developed a general framework in which the Schlesinger system and Korotkin and 
Samtleben's isomonodromic system are placed, along with generalizations to higher genus 
Riemann surfaces, in a unified way. Some more examples of matrix systems with different 
structures are also known jT21 EH HH 02] Compared with Okamoto and Iwasaki's for- 
mulation, these "elliptic analogues of the Schlesinger system" are obtained on an entirely 
different ground, such as conformal field theories, vector bundles on a torus, KZ equations, 
and (classical or quantum) integrable systems. This can be seen in the structure of the 
matrix linear system 



for which isomonodromic deformations are constructed. Namely, the matrix L(z) ("L- 
matrix" ) in these examples is borrowed from the isospectral Lax equation of an integrable 
system, though the Lax equation of isomonodromic deformations takes the non-isospectral 
form 



Each of those isomonodromic systems is thus accompanied by an isospectral partner. 

The correspondence between isospectral and isomonodromic systems will have a num- 
ber of significant implications. Among them, we are particularly interested in the role of 
"spectral Darboux coordinates". The notion of spectral Darboux coordinates was intro- 
duced by the Montreal group for isospectral systems with a rational L-matrix [16j . and 
extended to isomonodromic systems on a sphere [T7tll8j. As they demonstrated for those 
cases, one can construct a mapping from the Lax equation of this type of systems to a 
dynamical system of a finite number of points Pi, ... , Pn on the spectral curve 



though the spectral curve itself becomes dynamical in the case of isomonodromic de- 
formations. Spectral Darboux coordinates are the coordinates Ai, . . . , Ajv, fJ>i, ■ ■ ■ , Hn of 
the moving points Pk = (\k,fik)- These coordinates lead to "separation of variables" of 
isospectral systems. The most classical case is the so called Moser systems [TH]; separation 
of variable of those systems was worked out by Moser himself. Remarkably, as Harnad 
and Wisse pointed out [T%] . almost the same story repeats on the isomonodromic side, 





r = {(z,w) I det(W - L(z)) = 0} 
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except that separability is lost there. In particular, this shows an algebro-geometric in- 
terpretation of Okamoto's reformulation jE] of Garnier's work ^] on the 2x2 Schlesinger 
system. 

This paper presents a similar approach to one of the "elliptic analogues" of the 
Schlesinger systems, namely, the aforementioned isomonodromic system of Korotkin and 
Samtleben (in a slightly modified form). The isospectral partner of this isomonodromic 
system is the Calogero-Gaudin system [2H1I2I] f° r the SU(2) group. Separation of variables 
of the usual SU(2) Calogero-Gaudin system has been developed by Brzezihski [22J and 
Enriquez et al. [2H1 (including "quantum separation of variables" in the sense of Sklyanin 
[21]). Our method is more or less parallel to theirs, in particular, that of Brzezihski. Ac- 
tually, it is a rational (rather than elliptic) model of the SU(2) Calogero-Gaudin system 
that he considered. Thus we are to extend his method in two-fold ways — firstly, to an 
elliptic model (which is the subject of the work of Enriquez et al. as well), and secondly, 
to an isomonodromic system. 

A main outcome of our consideration (summarized in Theorem 1, Section V) is that 
the isomonodromic SU(2) Calogero-Moser system can be converted to a non-autonomous 
Hamiltonian system in the spectral Darboux coordinates. The Hamiltonians of this non- 
autonomous system turn out to be a considerably intricate functions of the Darboux 
coordinates and the time variables. Remarkably, however, a very similar Hamiltonian 
system has been discovered by Okamoto jH] for isomonodromic deformations of a second 
order scalar ODE on a torus. We shall show a natural explanation of this coincidence 
from our pouint of view. 

This paper is organized as follows. Sections II and III are for preparation. In Section 
II, the Poisson structure of the L-matrix of the Calogero-Gaudin systems is reviewed. In 
Section III, the isomonodromic system is formulated in terms of two canonically conjugate 
"Calogero variables" and a set of "spin variables". Section IV and V are the main part 
of this paper. In Section IV, the spectral curve and the spectral Darboux coordinates are 
introduced. In Section V, the non-autonomous Hamiltonian system is derived. Section 
VI deals with the relation to isomonodromic deformations of a second order scalar ODE. 
Section VII is for conclusion and supplementary remarks. Part of technical details are 
collected in Appendices. 
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II L-matrix and Poisson structure 



II. 1 L-matrix 

Following the idea of Korotkin and Samtleben JU|, we start from the L-matrix 

L(z) 



p * C{z-tj)A* faz-t^A] , (;|) 

\0 -P ) j^y^-^z-t^A} -((z-t^A* 

where q and p are Calogero variables, Af and Aj are spin variables, ((z) denotes the 
Weierstrass ( function and <f>(u, z) the auxiliary function that is widely used in the study 
of systems of the Calogero type: 

C(Z) = — ^y, (f>{u, Z) = ' 2 

Here a(z) is Weierstrass sigma function, and the prime stands for a derivative, i.e., <y'{z) = 
da(z)/dz. Let 2u± and 2uj 3 denote the primitive periods of the Weierstrass functions. 
Throughout this paper, we assume that tj ^ tk if j ^ k. 

This L-matrix is slightly different from that of Korotkin and Samtleben jl()j . They 
use Jacobi's elliptic theta function rather than Weierstrass' sigma function a. Their 
L-matrix is thereby more suited for formulating isomonodromic deformations against the 
modulus r. We dare to modify Korotkin and Samtleben's L-matrix because this simplifies 
the use of interpolation formulae of elliptic functions. It should be possible to start from 
the L-matrix of Korotkin and Samtleben and to derive substantially the same results, 
though we shall not pursue it in this paper. 

The Poisson structure of the dynamical variables is a standard one. The Calogero 
variables q,p are, in fact, the relative coordinate qi — q2 and momentum (pi — P2)/2 
of a two body system with canonical variables (qi,q2,Pi,P2), and become a canonically 
conjugate pair {q,p} = 1 in themselves. The spin variables Af, Aj obey the su(2) relations 

{A% A±} = ±5 3k At, {A+, A~ k } = 25 jh A* (3) 

with respect to the Poisson bracket. 

The Poisson bracket of the spin variables is nothing but the Kostant-Killilov bracket 
for the residue matrix 



/ 



A j 



A j A 3 



\ A 1 



(4) 



of L(z) at z — tj. The conjugacy class 



O j = {A j \A j ~diag(6 j /2,-e j /2)} (5) 

of semi-simple matrices with fixed eigenvalues ±0j/2 is a maximal (two-dimensional) 
symplectic leaf of this Poisson structure. One can use a canonically conjugate pair (xj, 
{xj,£j} = 1, to parametrize this symplectic leaf as follows: 



Note that this parametrization is consistent with the Poisson bracket of Af' 3 . 



II. 2 Poisson bracket of L-matrix elements 

Let us write the matrix elements of L(z) as 



L(z) 



A{u) B{u) " 



C{u) -A(u) 



More explicitly, 



where 



Uu, z) = = -<f>(u, z)(C(z -u) + C(«)). 



(7) 



N 

c(u) = 5>(-<z,*-ti)At. 

The non-zero Poisson brackets of these matrix elements take the form 

{A(z),B(w)} = B(z) ( f ) (-q,z-w)-B(w)((z-w), (8) 
{A(z),C(w)} = -C(z)d>(q,z-w) + C(w)az-w), (9) 

N 

{B(z),C(w)} = 2{A{z)-A{w))<P{q,z-w) + 2<p u {q,z-w)Y J A). (10) 

j'=i 



Thus the Poisson algebra of the matrix elements of L(z) almost closes up to the extra term 
proportional to X/JLi Ap which is later set to zero in order to derive the Lax equation. 

These Poisson commutation relations can be easily verified by direct calculations using 
the functional identity 



(u, z)(j)(—u, w) + 4>(u, z - w)(((z) - ((w)) + <j> u (u, z - w) = 



of the auxiliary functions. This functional identity is a consequence of the more general 
one 



(it, z)cj)(v, w) + 4>{u + v, z)(p(—v, z — w) — <fi{u + v, w)<p(u, z — w) = 0, 



;i2) 



from which the former identity can be derived by letting v — ► —it. 

The Poisson structure of the L(z)-matrix elements can be cast into the compact form 



{L(z) f L(w)} 



{L ab (z),L cd (w)}E ab ® E cd 



a,b,cA 



= [L(z) ® / + 1 ® L(w), r{z - w)} + 2 dr( ' Z W) ^ A% (13) 

where E ab denotes the matrix with the (a, b) element equal to 1 and the other elements 
vanishing. The r-matrix takes the form 

r(z-w) = ((z - w)E n (g) E u + <f>(q, z - w)E 12 ® E 21 

+ 4>(-q, z - w)E 21 (g) E 12 + C(z - w)E 22 ® E 22 
[ ((z-w) \ 

(j)(q,z-w) 

(f)(-q,z-w) 

V ((z-w) J 

which is a special case of the well known dynamical r-matrix of the elliptic Calogero-Moser 
system [23 123 HZ]- 



(14) 
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Ill Hamiltonians and Lax equations 



III.l Hamiltonians 

The fundamental Poisson commutation relation (|13|) implies that the standard quadratic 
Hamiltonians 

Hj = Res - Tr L(z) 2 (15) 
z—tj 2 

are not Poisson-commutative in themselves, but commute up to a term proportional to 

V iV A3. 

{n v n k }^j2 A3 j- ( 16 ) 

The factor J2jLi A? itself commutes with TC/s: 

N 

3=1 

To obtain a commuting set of flows, therefore, one has to impose the constraint 

JV 

E^ = 0. (18) 

3=1 

Note that J2?=i A) is an infinitesimal generator of the diagonal gauge transformations 

( e< . 

Aj i ► ^A^, p = . (19) 

I -e 7 

Because of the presence of the constraint (|T5jl . as Korotkin and Samtleben [TU] pointed 
out, one has to modify the naive Hamiltonians Tij by a term proportional to the left hand 
side of the constraint as 

N 

U 3 = Hj - 2 Pj J2 A\. (20) 

k=X 

The multipliers Pj, j = 1, . . . , N, are assumed to satisfy the relation 

N 

EPi = p> ( 21 ) 

3=1 
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which ensures the consistency of Lax and zero-curvature equations that we shall derive 
later on. A set of commuting flows are now defined on the reduced phase space by the 
canonical equations 

The Poisson brackets on the right hand side are understood to be calculated as 

N 

{F, H 3 } = {F, Hj} - 2p 3 {F, £ Aj}. (23) 

3=1 

Namely, we first calculate the Poisson bracket in the unconstrained variables, then impose 
the constraint. Since the term containing {F, pj} disappears upon imposing the constraint, 
we leave the Poisson brackets with p^s undefined. More explicitly, the equations of motion 
read as follows. 

1. Equations of motion of q and p: 

^- - 2A 3 

at j ~ 2A » 

= -Y,MQ,tj-t k )AfAT + J2M-Q,t J -t k )AjA+. (24) 

3 k^j k^j 

2. Equations of motion of A^: 

^ 13 A± , ni/i„ , + \ A± A 3 



T 2((t j ~ tk)A d At ± 2<j>(±q, tj - t k )AfAl ± 2 Pj A± (j + k), 



3 



dt 

f)A± 

-± = T 2 £ C(tk - U) Af A\ ±2j2<f>(Tq,t k -t i )AlAt T 2 J>A£. (25) 
otk t+k t+k i+k 



3. Equations of motion of A\: 

dA 3 

= -<f>(q, tj - t k )AjA^ + </>(-q, tj - t k )A~A% {j^k), 

f)A3 

-kt = E^t k -t e )AiAi -J2<P(-q,tk-t e )A k Al. (26) 

ol k t+k e^k 

In particular, the sum of these flows turns out to be trivial: 

™ 

j=l Ul 3 j=l Ul 3 j=l Ul 3 

The t-dependence of p/s cannot be determined in this way; we shall derive a set of 
differential equations for p/s later on. 
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III.2 Calculating {L{z),Hj} 

As an intermediate step towards Lax equations, we now consider the Poisson bracket 

{L(z),Hj} = {L(z),H j }-2 Pj {L(z),J2Al} ( 

k=i 

of L(z) with the modified Hamiltonians Hj. A clue is the general formula 

{L(z), -TrL(w) n } = Tr 2 ({L(z)®L(w)}I®L(w) n - 1 ), (29) 
where Tr 2 denotes the trace over the second component of the tensor product: 
Tr 2 ( X! X abcd E ab (g) E c dj = y^(y^X fflbec )-E a fe. 



a,b,c,d 



a,b c 



Plugging the Poisson commutation relation (fT3*j) into this formula, one obtains the identity 



{Hz), \ TrL(w) 2 | 



L(z), Tr 2 (r(^- w)I ® L(w))\ 
dr(z — w) x 



+ 2Tr 2 



dq 



-I®L{w) )Y,A). 



(30) 



3=1 



Extracting the residues at w — tj yields the Poisson bracket {L(z\ TCj}- Note that the 
second term on the right hand side disappears upon imposing the constraint. Since 



Res Tr 2 (r(^ — w)I <8> L{w)) 



w=t 



{-q,z-tj)Af -C(z-tj)A^ 



the Poisson bracket with Tij eventually takes the form 



where 



{L(z),Hj} — [L(z),Mj(z)], 



( 



Mj(z) 



{-q,z-t j )A+ -((z-t^A* 
On the other hand, the Poisson bracket with Y.k=i A\ can be readily calculated as 



(31) 



V 



A' 



N ( 

{L(z),J2Al} = £ 

k=l k=l 







z-t k )A k 
y<P{-q,z-t k )A+ 
[L(z),diag(-l/2,l/2)]. 
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One thus finds that 



{L(z),Hj} — [L(z),Mj(z)], (32) 

where 

Mj (z) = Mj (z) + diagfe , - Pj ) . (33) 

III. 3 Lax equations 

We are now ready to derive the Lax equations. 

By the Leibniz rule, the t-derivatives of the matrix elements of L(z) can be written as 

= E z - t k )^A- k + 0(g, z - t k )^) - 0'(g, z - tj )Aj, 

3 k=l 3 3 



and 



where 

^( U> = = _0( M , z)(C(M _ Z ) + C(z) ). 

Notice here that the terms containing (f)'(±q, z — tj) and ((z — tj) coincide with the matrix 
elements of —dMj(z)/dz. On the other hand, the 1,-derivatives of the dynamical variables 
q,p,Ak can be expressed as the Poisson bracket with Hj. Consequently, 

— -mz^Hj}-—^-. 

Combining this with the foregoing calculation of {L(z),Hj}, we eventually obtain the Lax 
equations 

of the isomonodromic type. Though we omit details, one can conversely derive the equa- 
tions of motion of q,p and A k s from these Lax equations. 
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As a remark, let us mention that these Lax equations are consistent with (|27|1 . This 
is a consequence of the linear relation (J2T|l among pj's and p. One can derive, from this 
linear relation, the linear relation 

N 

Y l M j (z) = L(z) (35) 
i=i 

among the matrices Mj(z) and £(z). The sum of the N Lax equations thereby yields the 
relation 

= "-&-■ (36) 

which is a restatement of (I2 - 



III. 4 Zero-curvature equations 

The Lax equations are self-consistent (i.e., define commutating flows) if and only if the 
zero-curvature equations 

dM k (z) dMj(z) 



()[ ()f + [M j (z),M h (z)]=0 (37) 

j k 

are satisfied. As we show below, these equations give a set of differential equations for 
Pi' 8 - 

The differential equations for pj's are obtained from the the diagonal part of the zero- 
curvature equations. We have only to consider the upper left component, because all ma- 
trices in the zero-curvature equations are trace-free. The diagonal part [Mj(z), M k (z)]n = 
— [Mj(z), Mk(z)]22 of the commutator can be rewritten as 

[Mj(z), Mk(z)] n = <P(q, z - tj)<j)(-q, z - t k )AjA+ - <j>(q, z - t k )<j>(-q, z - tj)A^A+ 
= -((j)(q,t k - tj)(C(z - tj) - C(z - t k )) + (J) u (q,t k - tj))AjA+ 
+ (0(g, t, - t k )(C(z - t k ) - C(z - tj)) + (j) u (q, tj - t k ))A k A+. 



The functional identity (fTTj) has been used in the last stage. The ^-dependent pieces are 
thus a linear combination of ((z — tj) and ((z — t k ). On the other hand, the derivative 
part dM ktU (z)/dtj — dMj^i (z) /dt k of the zero-curvature equation can be expressed as 

dM kM (z) dM^jz) dp k _ n M 
dtj dt k dtj + ^ {Z k> dtj dt k ^ h) dt k 
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Upon summation with the foregoing expression of the commutator part, the ^-dependent 
pieces ((z—tj)(- ■ -) + ((z — tk)(- ■ •) turn out to cancel out as a consequence of the equations 
of motion of Aj and A\. One is thus left with the differential equations 

§^ - ^ + MQ, tj - h)A+A- k - cj> u {q, t k - t^AjAt = (38) 

j k 

for pj's. 

One can examine the off-diagonal part of the zero-curvature equations in much the 
same way. This however leads to no new equations (see Appendix C). Namely, all equa- 
tions are satisfied under the equations of motion of q,p and Aj's. 

The final problem is the existence of a solution of ()38|) . To this end, it is convenient 
to convert ()38|) to an exterior differential equation. Let 9 denote the one-form 

N 

6 = Y,Pjdtj- (39) 

3=1 

f!38|) thereby turns into the exterior differential equation 

d6 = to, (40) 

where 

■y N y N 

j,k=l j,k=l 

N 

= E <t>(q,tj - t k )(C(tj -t h -q) + ((q))A+Aj;. (41) 
j,k=l 

Consequently, the Frobenius integrability of ([38)1 is equivalent to the closedness 

du = (42) 

of to. Actually, this integrability condition turns out to be satisfied under the equations 
of motion of q,p and Aj (see Appendix D). 

IV Spectral curve and Darboux coordinates 
IV. 1 Spectral curve 

The spectral curve is defined by the eigenvalue equation 

det{wl - L(z)) = 0. (43) 
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Since L(z) is trace-free, the left hand side can be rewritten as 

det(wI-L(z)) = w 2 + detL(z) = w 2 - ~ Tr L(z) 2 . (44) 

Under the constraint (|18|1. the matrix elements of L(z) enjoy the following quasi- 
periodicity along the period lattice of the torus: 

L{z + 2mcui + 2nu 3 ) = e - (mrn+nr *)Q L(z)e (m711+nV3)Q , (45) 

where Q is the diagonal matrix Q = diag(g, — q), and rji and r] 3 are the values of ((z) at 
z = uji,uj 3 . The quasi-periodicity of L(z) is a consequence of the quasi-periodicity of ((z) 
and (p(u, z), 

C(z + 2mu 1 + 2nuj 3 ) = ((z) + 2mr] 1 + 2noo 3 , (46) 
<f>(u, z + 2m Wl + 2mu 3 ) = e - 2mm ' 2nr,3 (f){u, z), (47) 

which are easy to confirm from the property of the sigma function. 

The quasi-periodicity of L(z), in particular, implies the double periodicity of Tr L(z) 2 /2, 
which thereby becomes an elliptic function with poles at z — t\, . . . , t^. Since 

= + O(l) 

Z Ij 

as z — > tj, this elliptic function has a double pole at z = tj with the leading coefficient 
equal to the quadratic Casimir 

C^\^A)J\ (48) 

of Aj. The residue is nothing but the Hamiltonian Tij. Thus Tr L(z) 2 /2 can be expressed 
as 

i N N 

- Tr L(z) 2 = £ C jP {z - t 3 ) + ^ HjCiz - tj) + Ho, (49) 

1 3=1 j=l 

where 7io is a constant term (which however depends on uj\ and uj 3 ). Also note that the 
Hamiltonians are not linearly independent, but obey the linear constraint 

N 

£Wi = 0. (50) 

3=1 
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This is a consequence of the the double periodicity of TrL(z) 2 /2. 

The spectral curve thus turns out to be a double covering of the torus. The branch 
points are located above the (possibly multiple) 2N zeros of TrL(z) 2 /2. If these zeros 
are all simple, the genus of the spectral curve is equal to N + 1. The spectral curve is 
time-dependent because of the extra term dMj(z)/dz on the right hand side of the Lax 
equation. 



IV. 2 Spectral Darboux coordinates 

The construction of spectral Darboux coordinates is parallel to the case of the rational 
(and isospectral) model |2~2"] . The "coordinate part" A 1; . . . , Ajv are defined as the N zeros 
(modulo the period lattice) of Li 2 (z), 

Li 2 {\ 3 ) = 0, (51) 
and the "momentum part" fix, . . . , //jv are defined to be the value of Ln(z) at these points, 

N 

to = MA,) =p + J2 C(Ai - t k )A\. (52) 

k=l 

In order to avoid a delicate problem, we assume throughout the following consideration 
that Xj 7^ Afc if j 7^ k. It is easy to see that (Xj, fij) sits on the spectral curve; L(Xj) takes 
the triangular form 

/ 

L(Xj) = 



to 
L2i(X j ) 




which implies that ±to are eigenvalues of L(Xj). 

The X/s are constrained by a linear relation. To see this, let us note that £12(2:) can 
be factorized as 



L l2{z) = K -±r, (53) 



where k is a constant that does not depend on z. The quasi-periodicity 
L 12 (z + 2mu 1 + 2nu 3 ) = e- {2mril+2nm)ci L 12 {z) 
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of Lviiz) implies that its zeros Ai, . . . , Xn are constrained as 



N N 

Xj — tj = q mod 2uiZ + 2CU3Z. 



(54) 



Since each Aj is defined only up to a difference by an element of the period lattice, let us 
redefine A/s, if necessary, such that this holds without "mod 2u>iZ + 2uj^Z": 



Of course this will be valid only for a local study of the system; this naive prescription 
has to be modified if one considers a global problem. 

Let us note here that the coefficient k transforms as k — > e 27 K under diagonal gauge 
transformation (|19|). One can thereby adjust k to any non-zero value, e.g., 



In other words, k is not a true dynamical degree of freedom. 
IV. 3 Time-dependent canonical transformation 

In order to prove the canonicity of these variables Xj, fij, we now restrict the spin variables 
onto the direct product 0% X ■ ■ ■ x On of the symplectic leaves and use the parametrization 
(|)J) by (xj,£j). Moreover, we tentatively relax the constraint (jl8j) . which now takes the 



N N 



(55) 



3=1 3=1 



K = I. 



(56) 



form 



N 



x & = °' 



(57) 



and restore it in the final stage. 

The factorization relation of L 12 (z) now reads 



If, ,,2 Q{Z) 

-\ i ^{q,z-t j )x j = K-^-y 



(58) 



where we have introduced the two functions 



N N 



Q(z) = l[a(z-X j ), P(z) = l[a(z-t j ). 



(59) 
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This reduces to the relations 

of the residues of both hand sides at z — tj. These relations show how the old variables 
Xj are connected with the new variables Xj (and tt). By logarithmic differentiation, these 
relations can be further converted to the linear relations 

N 

2dlogXj = d log k + d log ajtj — X k ) — d log ajtj — t k ) 

k=l k^j 
N 

= d\ogK+Y,t(tj-h)(dt j -d\ h )-Y l C{tj-tk)(dt j -dt k ). (61) 

k=l kjtj 

of differential forms. 



Our goal is to derive a relation between the canonical one-forms Y2j=i £,jdxj + pdq and 



Yuj=i HjdXj. We first multiply the both hand sides of the last relation by XjQ/2, sum over 
j = 1, . . . , N, and add pdq to both hand sides. We then obtain the linear relation 

N ^ N ^ N 

J2£jdxj+pdq = -^x j i j d\ogK + XjijC^tj- X^i^dtj- d\ k ) 
j=\ 1 j=l 1 j,k=l 

--^ -Os/UO - tk)(dtj - dt k ) + pdq. 
On the other hand, by differentiating (jo3j) . we have the relation 

N N 

dq = ^2 dXj — dtj, 

3=1 3=1 

which we can use to eliminate the differential dq on the right hand side of the foregoing 
linear relation of one-forms. The right hand side thereby becomes a linear combination 
of dlog/c, dX/s and dt/s, and the coefficient of dXj turns out to be equal to fij by (|52p. 
We thus eventually find that 

n y N N 

tjdxj +pdq = x£jd log k + ^3 d ^j 

3=1 3=1 3=1 

N ^ N 

-pE rft J + 9 x 3^c(tj - h)dtj 

3=1 Z 3,k=l 

-l^^Cfa-tkKdtj-dtk). (62) 

j¥=k 
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The last equation shows that Xj and fij are Darboux coordinates of the canonical 
one-form YljLi^jdXj + pdq, and that \ogn is a conjugate variable of the left hand side of 
constraint (|57jl . This interpretation is fully parallel to the spectral description of rational 
isospectral systems [TT) | IT7 | ITHj. 

An essential difference lies in the fact that the time variables explicitly enter the 
relation between the two canonical one-forms. This means that the spectral Darboux 
coordinates are connected with the old variables (xj,£j, q,p) by a time- dependent canonical 
transformation. Accordingly, the Hamiltonians (which is denoted by Hj in the following) 
in the spectral Darboux coordinates differ from the Hamiltonians TCj in the old variables 
( x j,£j,q,p). Their relation is to be determined by the fundamental formula 

N N N N 

ijdxj + pdq — £ Hjdtj = £ HjdXj — £ Hjdtj, (63) 
j=i j=i j=i j=i 

in which we have imposed the constraint (|57jl . More explicitly, the new Hamiltonians are 
defined as 

1 N 1 
Hj = H J +p- (J2 C(tj - A fc ) - E C(tj - t k )) +-E Cfe - t k )x k £ k . (64) 

fc=l kjtj kjtj 

The goal of the next section is to rewrite the right hand side in terms of the spectral 
Darboux coordinates. 

V Hamiltonian system in spectral Darboux coordi- 
nates 

V.l Linear equations characterizing Hamiltonians 

Let us recall that the pairs of the spectral Darboux coordinates all sit on the 

spectral curve. Therefore the equations 

N N 

A = E - tj) + £ HjCiXk - tj) + n 

j=l 3=1 

are satisfied for k = 1, . . . , N. These equations, along with the linear constraint (|5U|) . may 
be thought of as a system of linear equations that determine Tl/s. In fact, as we shall 
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discuss afterwards, these linear equations can be solved for TC/s, which thereby becomes 
an explicit function of the spectral Darboux coordinates (and of the time variables). 

If the system in consideration were an isospectral system (such as the Moser system or 
the usual Calogero-Gaudin system), the time variables would not appear here explicitly, 
and TC/s would be the Hamiltonians that we have sought. The only problem would have 
been to solve the foregoing linear equations for TLj. This is, actually, what Brzezihski 
[22] and Enriquez et al. (2H1 did in their work on separation of variables of the SU(2) 
Calogero-Gaudin system. 

In the present case, the true Hamiltonians are not Tt/s but H^s. We have to rewrite 
the extra terms on the right hand side of (|64|). too, as a function of the spectral Darboux 
coordinates. This is another problem that we have to solve. 

To this end, let us note that the defining equation (|52*|) of fij, which can be rewritten 

as 

N p 

^k = EC(A fc -ti)^+p, 

may be thought of as a system of linear equations for p and Xj^j/2. If one can solve these 
equations for p and Xj^/2, the solution should be an expression of p and Xj£j/2 as a 
function of the spectral Darboux coordinates. Remarkably, these linear equations have 
the same coefficients as the foregoing linear equations for H and TCj. Moreover, Xj^j/2 
obey the linear constraint (|57jl . in perfect analogy with the linear constraint (|5U|) for TCj. 

Thus the two problems, one for TCj and the other for the extra terms in (|64|). can be 
reduced to a single problem, namely, solving a system of linear equations of the form 

N 

Y / ah-t j )x j +x = b k (k = i,...,N), 

(65) 

i=i 

As we shall show below, this system of linear equations has a unique and explicit solution. 
V.2 Solution of linear equations 

We assume that q ^ mod 2a>iZ + 2u;3Z or, equivalently, a(q) ^ 0. The following ensures 
the uniqueness of solution of f)fc>5j) . 
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Lemma 1 If bj = ■ ■ ■ = = 0, then Xq = X\ = ■ ■ ■ = Xn = 0. 



Proof. Consider the function 

N 

f{z) = y £C(z-t j )X j +X . 

3=1 

The first N equations of (|fi5|) imply that f(z) has zeros at z = Xi, . . . , Xn- The remaining 
one ensures that f(z) is a doubly periodic meromorphic function on the z-plane. All 
possible poles are obviously simple and confined to z = tx, ■ ■ ■ , tjv an d their translations 
by the period lattice. Therefore, if f(z) is not identically zero, the zeros Xj and the poles 
tj are constrained as 

N N 

J2 X 3 ~Y1 t 3 = m ° d 2UJ 1 L + 2 ^3 Z ' 
3=1 3=1 

but this contradicts the assumption that q ^ 0; recall the constraint (|55j). Thus f(z) is 
identically zero, and all the coefficients Xq, X\, • ■ • , Xn have to be zero. Q.E.D. 

Having proven the uniqueness, the problem is to find a solution by any means. This 
can be done with the aid of an elliptic analogue of Lagrange's interpolation formula (see 
Appendix A). 

Lemma 2 A solution of h6fy) is given by 

x v Q^) P ^k)a{t 3 -X k + q)b k 

& P'(t j )Q'()*)*(t j - XkMq)' 1 ) 

x >P Qjt^pjXkMtj -x k + q )gx k - tj - q)b k 

k=i P'^Q'iX^-X^q) ■ 1 > 

Proof. We have only to confirm that these Xj and Xq do satisfy (jfi5j) . The last equation 
of (|65|) is indeed satisfied as (jA.l|) shows. As regards the other equations of (|65J1 . the main 
task is to calculate 

gC(A, - t,)X, - g(E p, fe)<Tfe _ Al) C(A, - W) Q^^-y 

We can use the two identities (|A.3j) and (|A.4|) to rewrite the sum over j on the right hand 
side, and find that 

£T ' J ' j ' " fc fc P%)Q'{X k )a{t 3 - X k )a{q) 

= -X + b h 
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which is nothing but the first N equations of (|65[). Q.E.D. 



V.3 Writing Hj explicitly 

Let us apply the foregoing formulae (pjj) and (joTjl of solution of (fBT>j) to the problem of 
deriving an explicit form of Hj as a function of the spectral Darboux coordinates. 
If we use the formulae to the case where 

N 

Xj = X = Ho, b k = fi 2 k — E Cip(Xk — ti), 

1=1 

we find the following expression of TCj and Tio: 

f Q(t j )P(X k )a(t j -\ k + q) , , f x 

^ " WWH^ ^-Egka, -;,)). (69) 

Similarly, if we use the formulae in the case where 

x j = 2 X o£oi x o =P, h = Hk- 

we find the following expression of and p as a function of the spectral Darboux 

coordinates: 

1 f ^ Q(t J )P(A fc )a(t J - X k + g)/^ 

2 ^ fe tl P'(t j )Q'(X k )(T(t j - X h )a(q) ' 1 J 

_ A Q(tj)P(X k )a(tj - X k + g)C(X k - tj - q)fi k 

p 3 r=i ntsWfrkMtj - x k )a( q ) ■ 1 } 

Thus we have been able to rewrite each term on the right hand side of (J64j) to an explicit 
function of the spectral Darboux coordinates. 

Although the extra terms on the right hand side of still appear to be in disorder, 
one can see by but straightforward calculations (see Appendix B) that the sum of these 
terms boils down to a form similar to the foregoing expression of Tif 

I N ^ i 

V ~ 2 X i^(E Cfe - X k) ~ E Cfe - t k )) + 7>J2 ((tj - tk)x k £ k 

k=l kj^j k^j 

h P'{t 3 )Q'{X k )a{t,-X k ) lUj q) Uj k,, ^ k 1 ; 
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Combining these results, we obtain the following expression of Hj in terms of the 
spectral Darboux coordinates: 

N 

x (l4 + (03 -h + q)~ 0j ~ A *)K - E c iP( x k ~ tij) . (73) 

i=i 

Eliminating q by (|55jl . we eventually obtain a final expression of the Hamiltonians: 

H y> g(t 3 -)P(A ife ) ( T(S^A f -E i#j tO 
J " jfe P'{tj)Q'(\k)v(tj - A fc )cr(EiIi A/ - tj) X 

A 7 " 

x (/4 + (Cfe - A fc + g) - Cfe - \k))Hk -£C,p(A fc - *»))■ ( 74 ) 

i=i 

In summary, we have proven the following. 

Theorem 1 The isomonodromic SU(2) Calogero-Gaudin system can be converted to the 
non-autonomous Hamiltonian system 



d\ k dHj dpik dHj 
dtj d/i k ' dtj dXj 

in the spectral Darboux coordinates Xj,fij. The Hamiltonians are given by J7^| ). 



(75) 



VI Relation to second order scalar ODE 
VI. 1 Deriving second order ODE 

The structure of the Hamiltonians Hj is very similar to Okamoto's Hamiltonians for 
isomonodromic deformations of a scalar ODE on a torus ||. This is not a coincidence, 
but can be explained in the same way as the case of the 2x2 Schlesinger system jH] . 
A clue is the fact that any 2x2 matrix system 

yields a second order scalar ODE of the form 

—r^ +Pl(z) — +P2{Z)yi = 0. (76) 
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The coefficients pi(z) and P2(z) of the latter are determined by the matrix L(z) as follows: 

Pl (z) = -TrL(z)-(logL 12 (z))', (77) 
p 2 (z) = detL(z)-L' u (z) + L 11 (z)(\ogL 12 (z)y. (78) 

In our case, L(z) is trace-free, so that the foregoing formulae of Pi(z) and p 2 (z) become 
slightly simpler: Firstly, P\[z) can be written explicitly as 

N N 

Pl (z) = -(logL 12 (z))' = - E C(z - Xj) + £ C(z - tj), (79) 

which implies that pi(z) is doubly periodic. Secondly, p-Az) is also doubly periodic (as 
the quasi-periodicity of the matrix elements of L(z) implies), and given by the formula 

p 2 (z) = -^TtL(z) 2 -L' n (z) + Ln(z)(\ogL 12 (z))'. 

One can see from this formula that p 2 (z) has simple poles at z — Xj and double poles at 
z = tj. Let us express p 2 (z) as 

N N N 

P2(z) = E «iP(« - tj) + E PiC(* - tj) + E 7iC(^ - Xj) + 5 (80) 
j=l j=l j=\ 

and determine the coefficients by Laurent expansion at the poles. 

1. The first coefficient ctj can be read off from the (z — tj)~ 2 term of the Laurent 
expansion of — Tr L(z) 2 /2: 

aj = -Cj = -TrA 2 /2. (81) 

2. The second coefficient f3j is the residue of p 2 (z) at z — tj. The term L' u (z) does not 
contribute to the residue. The residue of the other terms at z — tj can be expressed 
as 



1 

z=t~ 2 



Res-TrL(^) 2 = -Hj 



and 



ResLii(z)(logLi 2 (2;)) / = -p - - E ((tj ~ t k )x k ^ k 

1 kftj 



1 N 

+ (E c(tj - a,) - Eu'o - t k j). 



k=l k^j 
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As (JH1J) shows, the sum of these two quantities is equal to —Hj. Therefore 

(3j = fesp 2 (z) = -Hj. (82) 



3. The thrid coefficient 7, is the residue of p 2 (z) at z — Xj. Since 

— Res - TiL(z) 2 = 

z=\j 2 

and 

Res L n (z) (log L 12 (z))' = L n (\j) = /ij, 

Z — \j 

7j can be expressed as 

7j = Res p 2 (z) = Hj. (83) 

Z = \j 

Thus the Hamiltonians Hj and the "momenta" fij can be identified with the residues 
of p 2 (z). Exactly the same relation can be seen in the case of Gamier 's isomonodromic 
system on a sphere [U]. 

VI. 2 Another form of second order ODE 

Strictly speaking, however, the second order ODE above differs from that of Okamoto [8] 
and Kawai j^j, who consider a linear ODE of the form 

^+p(z)y = 0. (84) 

At least formally, this discrepancy can be removed by the "gauge transformation" 

yi = exp(-iy" Pl ( z )dz)y. (85) 

The coefficient p(z) is given by 

P(z) = -\v'Az) ~ \pi{z) 2 +p 2 {z). (86) 

Note, however, that this is actually a delicate procedure, because the gauge transformation 
might spoil the isomonodromic property. Fortunately, the present case is free from this 
problem: The gauge transformation takes the form 
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and since the factor in front of y has constant monodromy, the isomonodromic property 
is preserved by the gauge transformation. 

The zero-th order term p(z) of the transformed ODE is a doubly periodic meromorphic 
function with second order poles at z — Xj,tj. The residues of p(z) at these poles can be 
readily determined: 

JV 

Resp(z) = /i,-2^C(A J -4) + 2^C(A, -A fc ), (88) 
z ~ 1 k=i k^j 



N 

Resp(z) = -H j -2'£C(t j -Xk) + '2^C{tj-tk). (89) 
■ ' k=i k+j 



It is rather these quantities that Okamoto [S] and Kawai jH] use as Hamiltonians and 
conjugate variables of Ay's. We can indeed reformulate our Hamiltonian system in that 
way. Namely, if we define 

Uj = Res p(z), Kj = —Resp(z), (90) 
these quantities satisfy the equation 

N N N N 

HjdXj — V] Hjdtj = ^ VjdXj — V] Kjdtj + exact form. (91) 

j=l 3=1 3=1 3=1 

This implies that Aj's and v/s are a new set of Darboux coordinates, and that the previous 
Hamiltonian system is now converted to the new Hamiltonian system 



d\ k dKj dv k dKj 
dtj dv k ' dtj dXj 

VI. 3 Reconstructing 2x2 matrix system 



(92) 



Let us now consider the inverse problem. Namely, given the isomonodromic deformations 
of the second oder scalar ODE above, the problem is to reconstruct a 2 x 2 matrix system. 
A similar problem is discussed by Okamoto [HI Section 3] in the case of isomonodromic 
deformations on a sphere. In our case, the presence of diagonal gauge transformations 
(|19|) allows us to fix the coefficient 1. Therefore we have only to show how 

to reconstruct the dynamical variables q,p,A k s of the matrix system from Xj,fij of the 
scalar ODE. 

One can indeed reconstruct the L-matrix L(z) by an algebraic procedure as follows. 
The first step is to reconsider (JHHJ) and (JHSJ as definition of Li 2 (z) and q. The coefficient k 
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is chosen to be k = 1 as remarked above. Aj's are thus determined. Secondly, let Ln(z) 
be a function of the form 

N 

L 11 {z)=p + J2<t>{Q,z-t j )A 3 j (93) 
i=i 

that satisfy the interpolation relations Lu(Xj) = fij for j = 1, . . . , N. As we have seen in 
Section V, these relations can be solved for A? and p under the constraint J2?=i ^| = 0. 
As (0) suggests, A+'s are to be determined as 



= gi^g (94) 



One can thus reconstruct L(z) from the scalar ODE. 



VII Conclusion 

We have applied the method of spectral Darboux coordinates to Korotkin and Samtleben's 
isomonodromic system on a torus ^U]. The isomonodromic system has been thus con- 
verted to a non-autonomous Hamiltonian system in the spectral Darboux coordinates. 
Although the Hamiltonians turn out to be a considerably intricate function, the method 
we have used is a rather straightforward analogue of the usual method for isomonodromic 
deformations on a sphere. 

Our non-autonomous Hamiltonian system may be thought of as an elliptic analogue of 
Garnier's isomonodromic systems EJ E] ■ Almost the same system has been derived by 
Okamoto from isomonodromic deformations of a second order scalar ODE on a torus jS] . 
We have seen how these two systems are related. Speaking differently, our approach from 
a 2 x 2 matrix system reveals a hidden algebro-geometric meaning of the Hamiltonian 
structure in Okamoto's work [UJ. 

An important lesson of the present work is that the notions of spectral curve and 
spectral Darboux coordinates persist to be useful and essential beyond isospectral defor- 
mations. This observation lies in the heart of the work of Harnad and Wisse [THj . We 
have confirmed it for an example of isomonodromic deformations on a torus. 

In this respect, an interesting problem is to describe the isomonodromic SU(2) pure 
Gaudin system H3] from the same point of view. Separation of variables of the 
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isospectral partner has been studied by Sklyanin and Takebe j28j (see also the paper of 
Hurtubise and Kjiri [2H| for geometric aspects). The work of Sklyanin and Takebe shows 
that separation of variables of this system is technically far more complicated than the 
Calogero-Gaudin system. This will be also the case for the isomonodromic analogue. 

Let us conclude the present consideration with a remark on trigonometric and ratio- 
nal analogues. The trigonometric and rational analogues of Korotkin and Samtleben's 
isomonodromic deformations can be obtained by replacing the basic functions cr(z), ((z) 
and 4>{u, z) by the following trigonometric or rational functions: 

1. Trigonometric model 

. . . . , . COS Z ,, . COS Z COSM . 

a( z) — sm z, Gz) = , (p{u J z) = . (95) 

sin z sm z sm u 

2. Rational model 

a(z) = z, C(z) = -, <j>(u,z) = ---. (96) 
z z u 

A hyperbolic model will be obtained if one replaces the trigonometric functions by the 
corresponding hyperbolic functions. These are nothing but the well known pattern of 
degeneration of the Calogero-Moser systems; the Calogero-Gaudin systems, too, obey 
this pattern. In fact, it is the rational model in this list that Brzezihski considered in his 
work [22] • One can formulate an isomonodromic partner of these degenerate Calogero- 
Gaudin systems as in the case of the elliptic model. Presumably, those isomonodromic 
systems will not be known in the literature. 
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A Interpolation formula 

Let us examine the auxiliary function 

Q{z)a(z - X k + q) 



h{z) 



P(z)a(z - A fc ) 
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This is a doubly periodic meromorphic function with simple zeros at Xj (j ^ k) and Xk — q 
and simple poles at t j (j — 1, . . . , N). By the residue theorem, the residues 

at the poles z = tj obey the sum-to-zero constraint 

Q{tj)a{tj -Xk + q) = n rAn 
P%)a{t 3 - X k ) 1 • > 

Let us consider the linear combination 

" Q( tj )a(tj-X k + q) 
£ P%)a{t 3 -X k ) Q[Z tj) 

of ((z — tj) weighted by these residues. Since this function is a doubly periodic meromor- 
phic function with the same set of simple poles and residues as f k (z), it differs from f k (z) 
by at most a constant: 

Q{z) °[ z 7 \\ q) = f Q^Mh-Xk + q) _ congtant 

P(z)a(z-X k ) £ P'(tj)a(tj-X k ) ^ jJ 

Moreover, since the left hand side vanishes at z = X k — q, the constant term on the right 
hand side can be easily determined as follows: 

constant = -g p>(t Mh -x k) c(Afc " *' - q) - 

We thus obtain the interpolation formula 

Q(z)a(z -X k + q) * Q(tj)a(t 3 - X k + q) 
P(z)a{z - X k ) j=l P\tj)a{tj - X k ) 

One can derive the following three identities from this interpolation formula. 

1. Since the left hand side of the interpolation formula vanishes at z — A; (I ^ k), 

£ Q pufX ^xf (C(Ai - tj) - C(A fc - t 3 - - g)) = (Z^fc). (A.3) 

2. By letting z — > A& in the interpolation formula, 

f gftx^-Afc + g) _ gWg) 
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3. By replacing k I, j —> k and separating a term from the sum, the interpolation 
formula takes the form 

\ - Q{t k )a(t k - \i + q) ^ _ 

Q(z)a(z -Xi + g) Qjt^ajtj - X t + g) _ 

P(z)a(z-X l ) P'itMtj - X t ) U " j) 

By letting z — > tj, 

\- Q(t k )o(t k - Xi + q) _ 

+ E^fe^^C(A,-4-g). (A.5) 
£r x r '{t k )a{t k - Xi) 



B Calculation of extra terms in (1641) 



Let us use (j70|) to rewrite the last piece on the right hand side of as 

2 Cfe ~ tfc)Xfc6 " % P'{t k )a{t k - Xi) ^ - 4) J QW^)- 
The sum over k ^ j arising here has been partially calculated in (jA.5|) . Using the identities 



on the right hand side of ()A.5jl . one can rewrite the foregoing quantity as 
1 N 

9 E C{tj - tk)Xh(k 

z k=l 

^ Q(t J )P(X l )a(t k -X l + q) t i e f _ , u f >, x , 

fr[ P\tj)Q'{Xi)a{tj - Xi)a{q) 

Q(t k )P(Xl)(T(t k ~Xj+ q)C(Xl - tfe - gVj 

^ P'(t k )Q'{Xi)(T(t h - Xi)a(q) 
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By (|70p and ()7ip. the first and third lines on the right hand side turn into the following 
form: 

N 



first line = --X& C(*j ~ *fc) ~ C(*j ~ V 

fc^j fc=i 



third line = —p. 

Since the sum of these two cancels the second and third pieces on the right hand side of 
(|51jl. we eventually obtain the identity 

1 N 1 
P - 2 X &{12 C(tj - A fe ) - Cfe - **)) + g £ Cfe - 
fe=i fc^j fc^j 

" QMPMafa — Xi + q) 
~ g PW(A^ - A.)o-(g) (Cfe - A ' + q) ~ ^ ~ A ' )K 

which is nothing but ()72|1. 



C Zero-curvature equations in more detail 

Since the diagonal part of the zero-curvature equation (J37)l has been specified in Section 
III, let us now examine the off-diagonal part. 

As regards the upper right part, the matrix element of the commutator [Mj(z), M k (z)\ 
reads 

[Mj(z),M k (z)] 12 = 2( Pj + az-t J )A 3 1 ) ( f ) (q,z-t k )A^ 

- 20(g, z - t 3 )A~( Vk + ((z - t k )Al). 

One can use the functional identity in the form 

cj)(u, z — w)((z) = — 4>{u, z)(f)(u, —w) + (f>(u, z — w)((w) — (fi u (u, z — w) 
to eliminate ((z — tj)<p(q, z — t k ) and <p(q, z — tj)((z — t k ) as 

[Mj(z), M k {z)} 12 = 24>(q, z - t k ){p 5 Al + ((t k - tj)J§Al + 4>(q, t k - t,)AjAl) 

- 2(f>(q, z - tj)(p k Aj + C(tj - t k )AjAl + <j)(q, tj - t k )A)Al) 

- 2<f) u {q, z - t k )A)Al + 2<f) u (q, z - tj)AjA\. 
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On the other hand, the derivative part of the zero-curvature equation can be expressed 
as 

dM k>12 {z) 8M hl2 {z) 
dtj dt k 

, . . dAl , , s dq ._ .. . dAj dq . _ 

3 3 k k 

Thus the (j) u terms cancel out in the zero-curvature equation and one is left with an 
equation of the form 

<t>{q, z - tj){- ■ •) + <t>{q,z - t k ){- ■ •) = 

with the coefficients (• ■ ■) that do not depend on z. As it turns out, these coefficients are 
exactly the same as some of the equations of motion of q, Aj and A k . One can thus see 
that this part of the zero- curvature equations is automatically satisfied. 

The lower left part of the zero-curvature equation can be treated in the same way. 



D Integrability of (SSj) 

In components, the closedness condition du = reads 

8 tj (4>(q, t 3 - t k ) (((t, -t k -q)+ C(q))A+A^) 

+ (cyclic permutations of j, k, £) = 0. (D.6) 

The goal is to show that these equations are indeed satisfied under the equations of motion 
(1231) . (f25jl. (127)1) of q,p and A's. Applying the Leibniz rule to the left hand side of (|TX6|l 
yields such terms as 

-<P(q,t 3 - * fc )((C(*i ~tk-q) + C(q)) 2 + C(tj ~ h - q) + C (q))^-A+ A k 

f)A+ r) A~ 

+<P(q,t 3 - t fc )(C(*i ~t k -q) + C(q))(-^-Aj; + 

and their cyclic permutations. One can eliminate the derivatives of q,p and A's by the 
equations of motion; this in turn yields linear and quadratic combinations of 0's. As 
regards the quadratic combinations, one can use the functional identity (fTT|) to reduce 
such terms to a linear combinations of 0's, e.g., 

<f>(q,tj - t h )<f>(q,ti - tj) = <f>(q,t e - t k )(C(q - U + t k ) - ((q) + ((tj - t k ) + ((t £ - tj)}, 
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etc. After some more algebra, one can thus eventually convert (jD.6|) to equations of the 
form 



4>(q, tj — tk)cjkeA~j~A k A\ + (cyclic permutations of j, k, £) — (D.7) 

where 

Cjkt = (C(t j -t k -q) + ({q)) 2 + C{t j -t k -q)+('(q) 
+ (C(tj ~t k -q) + ((q))(((te - t,) - C(t t - t k )) 

+ (C(t £ -t k -q) + C(g))(-C(*i -tk-q)- C(q) + C(k - h) + ((t, - u)) 

+ (Cfe ~U-q) + C{q)){-C{tj -t k -q)- C(q) + Cfe " U) + ((U - t k )). 

Actually, these coefficients Cj k £ turn out to vanish identically. One can indeed verify that 
the function 

f jM (z) = (((t j -t k -z) + ((z)) 2 + C(t j -t k -z) + ('(z) 

+ (C(tj -t k -z) + ((z))(((u - tj) - C{U - t k )) 

+ {C{te -t k -z)+ C(z))(-C(tj -t k -z)- ((z) + C(U - t k ) + Citj ~ U)) 

+ (C(tj -U-z) + C(z))(-C(tj -tk-z)- C(z) + ((tj - t t ) + c(t t - t k )) 

of the complex variable z is a doubly periodic entire function with a zero at z = (tj—t k )/2; 
this implies that fj k e(z) = 0. 
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